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We performed the combine fits of SLAC, BCDMS and NMC data at NLO level. 
The data are fitted very well with different parametrizations of higher twist con- 
tributions. The theorety-based parametrizations lead to the low values of QCD 
\S~ A ~, 1 /in 4- 01 WoT/ f™. fn/f2> 
MS 



parameter Ai_Z ~ 140 ± 21 MeV (or a s {M%) ~ 0.103 ± 0.002). 



The accuracy of the present data for deep inelastic (DIS) structure func- 
tions (SF) reached the level at which the <5 2 -dependence of logarithmic QCD- 
motivated and power-like ones may be studied separately (see, for example, 
the recent reviews u and references therein) . In the present letter we analyse 
at NLO order of perturbative QCD the mpst known DIS SF ^2(2;, Q 2 ) taking 
into account SLAC, BCDMS and NMC □ experimental data. We stress the 
power-like effects, so-called twist-4 (i.e. ~ l/Q 2 ) and twist-6 (i.e. ~ 1/Q 4 ) 
contributions. To our purposes we represent the SF F 2 (x,Q 2 ) as the contri- 
bution of the part F^ 00 (x,Q 2 ) (including target mass corrections (TMC)) 
described by perturbative QCD and the nonperturabative part: 



Contrary to standard fits (see, for example, i) when the direct numerical 
calculations based on DGLAP equation are used to evaluate of SF, we use the 
exact solution of DGLAP equation for the Mellin moments of (x, Q 2 ) 

and the subsequent reproduction of F 2 (x n Q 2 ^ at some Q 2 -value with help of 
the Jacobi Polynomial expansion methodmjtj 

F 2 N —{x 1 Q 2 )=x a (\-xf £ e^(x)5>5 n W)A^ +2 (Q 2 )> ( 2 ) 

n=0 j=0 

where Q^'^ are the Jacobi polynomials and a, (3 are their parameters, fixed 
by the condition of the requirement of the minimization of the error of the 
reconstruction of the structure functions. 
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For n-space eq. (|T|) transforms to 



■M™(Q% (3) 

where (hereafter (a = 4,6)) C a (n) = x n - 2 h a {x)dx/M^ s {Q 2 ). The Q 2 - 
evolution of the moments M^ S (Q 2 ) is given by the well known perturbation 
QCD formulae (seeQ). The moments MP^ CD (Q 2 ) are distinct from M^ S (Q 2 ) 
by including TMC. The moment M^ S (Q^) is theoretical input of our analy- 
sis and its parameters should be found together with h±(x), h%{x) and QCD 
parameter A by the fits of experimental data. 

The shapes h a (x) of the higher twist (HT) corrections are of primary 
consideration in our analysis. They are choosen in the several different forms. 

• The twist-4 and twist-d terms are fixed in agreement with the infrared 
renormalon model (seeQ and references therein \_ The values of HT coef- 
ficient functions Ci{ri) and C&{n) are given (seeH) in the form C^ en (n) = 
B a (n) ■ A^ en , where B a {n) are known n-dependent functions and Aa en 
are some coefficients. The values of A^ en have been estimated (inu) as 
A r f n = 0.2 GeV 2 and A r 6 en = {A\ en ) 2 = 0.04 GeV 4 . We use them here 
as free parameters. 

• The twist-4 term in the form h 4 (x) ~ -^lnF 2 NS {x, Q 2 ) ~ 1/(1 - x). 
This behaviour matches the fact that higher twist effects are usually 
important only at higher x. The twist-4 coefficient function has the form 
Ct r (n) = (n-l)Af er . 

• HT terms are considered as free parameters at each x\ bin. They have 
the form hl ree (x) = X)i=i ha{xi), where / is the number of bins. The 
constants h a (xi) (one per a:- bin) parametrize x-dependence of h' a ree {x). 

To clear up the importance of HT terms we fit SLAC, BCDMS and NMC 
(H 2 and D 2 ) dataEl (including the systematical errors), keeping identical form 
of perturbative part at NLO approximation. We choose hereafter Q 2 = 10 
GeV 2 , N max — 6, 0.3 < x < 0.9 and 1 = 6. Basic characteristics of the quality 
of the fits are x 2 /DOF for SF F 2 and for its slope dlnF 2 j 'dlnQ 2 , which has 
very sensitive perturbative properties (seeQ). We use MINUIT program for 
minimization of two \ 2 values: 

2 ( h = dlnF A 

' ^ dlnQ 2 ) 

We obtain the following results: 
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M n (Q 2 ) = M^ CD (Q 2 



C 4 (n) , C 6 (n) 



Q 2 



Q 4 



x 2 (F 2 ) 



jpteor 



Ap exp 



and x 2 (slope) 



foexp ^te 

Ab ex P 



• Only perturbative QCD part is included (i.e. C a = 0, M 2 ucl = 0): 

qo^S 1 Q77 

x 2 /DOF(F 2 ) = — , X 2 /£>OF(slope) = Ag = 300 ± 5 MeV 

• Only TMC are included (i.e. C a = 0): 

X 2 /DOF{F 2 ) = x 2 /£OF(slope) = Ag = 257 ± 5 MeV 

• TMC and h{ ree (see Table) are included: 

X 2 /DOF(F 2 ) = |||, X 2 ADOF(slope) = ^ A gJ = 287 ± 28 MeF 

• TMC, Cf r and ^ ree (see Table) are included: 

X 2 /DOF(F 2 ) = ||, x 2 / J DOF(slo P e) = ™, Ag = 154±22M e y 
A der (H 2 ) = 0.244 ± 0.041 Ge^ 2 , v4 rfer (D 2 ) = 0.162 ± 0.037 GeF 2 

• TMC, C£ e ™ and C\ en are included: 

560 1 5 

X 2 /DOF{F 2 ) = — , X 2 /^OF(slope) = y , Ag = 125 ± 17 MeV 

A r 4 en (H 2 ) = 0.485 ± 0.035 GeV 2 , A r 4 en (D 2 ) = 0.457 ± 0.032 GeV 2 
A r 6 en (H 2 ) = 0.0070 ± 0.0005 GeV 4 , A r 6 en (D 2 ) = 0.0024 ± 0.0006 GeV 4 

Let us now present the main results. 



The different parametrizations oLHT terms lead to very good fits of 
SLAC, BCDMS and NMC data H. Including twist-4 term correction 
as free parameters at each xi bin leads to A£l^ ~ 286 ± 28 MeV (or 
a s (M|) ~ 0.114 ±0.002). However, the HT terns, which are obtained in 
the framework of infrared renormalon model and also in the form of the 
derivative 4-lnF^ s {x, Q 2 ), lead to the lower values of QCD parameter 
Ag - 140 ± 21 MeV (or a s (M§) ~ 0.103 ± 0.002). It is in perfect 
agreement with the recent results U of DELPHI Collaboration, where the 
strong reduction of the a s (M|) value was observed when renormalon- 
type power corrections have been included. The reduction in □ is even 
stronger becaus e th e nonperturbative part for the cross sections of e + e~ 
collisions starts E3 with the l/^/Q 2 term. Moreover, similar effect of 
the reduction of the a s (M^) value has been also observed earlier in the 
analyses of DIS data based on NNLO approximation Q af perturbative 
part and also on scheme-invariant perturbation theory □ resumming a 
part of higher order terms. 
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• Twist-4 terms are very similar to ones obtained ini, when we considered 
them as free parameters. 

• Twist-6 terms are different for proton and deutron SF that seems to 
indicate a nonzero nuclear correlations in deutron. 



Table. The values of the parameters for h{ ree ( GeV 2 ) and h^ ree 



GeV 4 





0.35 


0.45 


0.55 


0.65 


0.75 


0.85 


h{ ree (H 2 ) 


-0.245 ± 
0.025 


-0.171 ± 
0.045 


-0.029 ± 
0.076 


0.336 ± 
0.133 


0.759 ± 
0.239 


1.41 ± 
0.60 


h{ ree (D 2 ) 


-0.205 ± 
0.030 


-0.064 ± 
0.048 


0.106 ± 
0.078 


0.271 ± 
0.128 


0.426 ± 
0.213 


1.44 ± 
0.55 


hi ree (H 2 ) 


-0.266 ± 
0.018 


-0.425 ± 
0.054 


-0.618 ± 
0.134 


-0.448 ± 
0.414 


-1.61 ± 
0.99 


-9.52 ± 
4.32 




-0.198 ± 
0.020 


-0.139 ± 
0.057 


0.003 ± 
0.137 


0.147 ± 
0.398 


-0.651 ± 
0.821 


3.72 ± 
4.36 
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